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Abstract 

We consider formal quantum Hamiltonian of a charged particle on the Poincare 
, disk in the presence of an Aharonov-Bohm magnetic vortex and a uniform magnetic 

I field. It is shown that this Hamiltonian admits a four-parameter family of self- 

adjoint extensions. Its resolvent and the density of states are calculated for natural 
values of the extension parameters. 
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1 Introduction 



CN ■ Quantum dynamics on the Poincare disk has long been a subject of theoretical inter- 
. est, mainly because of the insights its study provides into the theory of quantum chaos. 
\0 ', Analyzed examples include, for instance, the free motion under the action of constant 
^ I magnetic fields [H [9l [23] , the Kepler problem [25] , the scattering by the Aharonov-Bohm 
O ' (AB) [211 [26] and Aharonov-Bohm-Coulomb [29] potentials, the study of point interac- 

tions [31 [5] and quantum Hall effect [6]. 

In the present paper, we consider the Hamiltonian of a charged spinless particle moving 
P_i! on the hyperbolic disk, pierced by an AB flux, in the presence of a uniform magnetic 
field. First part of this work is rather standard: we determine the admissible boundary 
conditions on the wave functions, using Krein's theory of self-adjoint extensions (SAEs) 
g ■ [I]- It turns out that in the most general case the formal Hamiltonian has deficiency 

indices (2, 2) and thus admits a four-parameter family of SAEs. Let us remark that 
similar results on the plane have been found in [21 E] in the case of zero magnetic field, 
X I and in [21] for non-zero fields; SAEs of the Dirac Hamiltonian on the plane have been 
c5 ; studied in [T6l[22]. 

The rest of this paper is devoted to the study of a particular extension, corresponding 
to the choice of regular boundary conditions at the position of the AB flux. We start by 
constructing certain integral representations for common eigenstates of this Hamiltonian 
and the angular momentum operator. These representations then allow to sum up the 
contributions coming from different angular momenta to the resolvent kernel, and to 
evaluate this kernel and the density of states in a closed form. 

The material is organized as follows. In Section 2, we introduce basic notations and 
study elementary solutions of the radial Schrodinger equation on the Poincare disk. Self- 
adjointness of the full AB Hamiltonian is discussed in Section 3. In Section 4 we find 
a compact expression for the resolvent of the regular extension (formulas I4.14p . fl4.19p - 
fl4.22p ). These relations represent the main result of the present work. The density of 
states, induced by the AB flux in the whole hyperbolic space (see f l5.8p - fl5.10p ). is obtained 
in Section 5. Some technical results are relegated to the appendices. 



1 



2 Free Hamiltonian on the Poincare disk 



2.1 Basic formulas 

Let us identify the Poincare disk D = SU{1, 1)/S0{2) with the interior of the unit circle 
|zp < 1 in the complex plane, equipped with the metric 

ds =gzzdzaz = R ^ (2.1) 

(1 - |Z|2) 

of constant Gaussian curvature — We consider a spinless particle moving on the 
disk and interacting with a magnetic field. The latter can be introduced as a connection 
1-form 

A = Azdz + As dz 

on the trivial t/(l)-bundle over D. Quantum dynamics of a particle of unit charge is 
described by the Hamiltonian 

H = -— {D,,D,}, (2.2) 

where = dz+ iAz and = 9^ + iAz are the usual covariant derivatives. To unburden 
formulas, we put the particle mass equal to 1/2 and h = c = 1 throughout the paper. 
In the remainder of the present section, the following vector potential is considered: 

iBR'^zdz — zdz 

It generates a curvature 2-form J-'^^\ proportional to the invariant volume measure djj = 
- gzz dz A dz. Indeed, we have 

=dA^''^ = Bdfi. 

Therefore, the potential (12.31) describes a uniform magnetic field of intensity B. Introduc- 
ing polar coordinates z = re*'^, z = re"**^, one can write the corresponding Hamiltonian 
as 

H^^^ = ld.r+'-dr + \d,, + ^^d,- -f^ A . (2.4) 

[ r 1 — 4 (1 — r^) J 

Note that the domain of H^^^ is not yet specified. It will be fixed in the next section by 
the requirement for the Hamiltonian to be a self-adjoint operator. According to Stone's 
theorem, this condition ensures the existence of consistent dynamics. 
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2.2 Radial Hamiltonians 

Formal Hamiltonian H^^^ commutes with the angular momentum operator L = —id^. 
Therefore, it leaves invariant the eigenspaces of L, spanned by the functions wi{r)e^^'^ 
(/ e Z). Being restricted to the eigenspace of L, characterized by the angular momentum 
Z, the Hamiltonian acts as follows: 

wi{r) ^ Hiwi{r), 

(2.5) 

Here we have introduced instead of -B a dimensionless parameter b = BR'^/4. 

It will be useful for us to let the parameter / to take on not only integer, but also 
arbitrary real values, and to study in some detail the properties of solutions of the radial 
Schrodinger equation 

'Hi - e) wi = 0. (2.6) 



In what follows it will be always assumed that k'^ G CyM"*" U {0}. It is also convenient to 
introduce instead of r a new variable t = r"^. 

We are interested in the solutions of (12. 6p leading to square integrable (with the 
measure ci/i) functions on D. These solutions should be then square integrable on the 

open interval / = (0, 1) with the measure d^t = ^^^^ — For each / G M there exists 

only one solution of (12. 6p . which is square integrable in the neighbourhood of the point 
t = 1. Its explicit form is 

wf\t) = t-'/^l-tr,F,ix-b,x + b-l,2x,l-t)= (2.7) 

where 



_ 1 + VI + 462 - A;2i?2 
X - ^ 

and 2Fi{a, z) denotes Gauss hypergeometric function. The branches of square roots 
are defined so that they take on real positive values for purely imaginary k. 

Similarly, for each / G (— cx), — 1] U [l,C)o) there is only one solution of (12. 6p . which 
is square integrable with respect to dfit near the point t = 0. The form of this solution 
depends on whether / > 1 or Z < —1. In the first case, i. e. for Z > 1, it is given by 

^(^^•+) (t) = tV2 (1 _ t)x + ^ _ ^ + 1 + (2.8) 

while for / < — 1 this solution is written as follows: 

wi'''~\t) = t-'l^ (1 - ,F,{x-b,X + b-l,l-l,t). (2.9) 

Note that for |/| < 1 both functions wl^^'^\t) are square integrable in the vicinity of the 
point t = and solve the radial Schrodinger equation (12. 6p . These solutions are linearly 
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independent except for / = 0. However, in the latter case the equation fl2.6p still admits 
two distinct solutions that are square integrable as t ^ 0: 

where u and v are any two linearly independent solutions of the hyper geometric equation 
with parameters a = x + b, l3 = x~b,'~f = l (one can choose them, for instance, according 
to the formulas 15.5.16 and 15.5.17 of [Ij). 

Let us now show that the solutions Wi^\t) and Wi^^'~^\t) are linearly independent for 
/ > — 1, and the solutions wl^\t) and w^^^' \t) are linearly independent for / < 1. This 
can be done by an explicit computation of their Wronskian 

W{h,f2) = fi-dtf2-dj,- h. 

Namely, using the connection and analytic continuation formulas for hypergeometric func- 
tions [1], one obtains 

w{wP{t),wl'''^\t)) = {tC^^,)-\ (2.10) 

with 

, r(x±WxT^. (2.11) 



Therefore, for k"^ G C\M"'" U {0} and |/| > 1 the equation fl2.6p has no square integrable 
solutions (with the measure dfit) on the whole interval /. This is true, in particular, for 
all radial Hamiltonians Hifzz of the free particle in a uniform magnetic field, except for 
the s-wave Hamiltonian Hq. In the case |/| < 1 the equation (12. 6p has exactly one square 
integrable solution, given by the formula (12. 7p . 

Let us now restrict the domain of Hi to 'D{Hi) = C^{I), i. e. to smooth compactly 
supported functions. Then the above remarks imply that 

• Hi is essentially self-adjoint for |/| > 1, 

• for |/| < 1 the operator Hi has deficiency indices (1, 1) and thus admits a one- 
parameter family of self-adjoint extensions (SAEs). 

Different extensions Hi'^^ (|/| < 1) are in one-to-one correspondence with the isometries 

between the deficiency subspaces JCf" = ker (^Hi =f ie^ , where e G may be chosen 

arbitrarily. They can be labeled by a real parameter 7 G [0, 27r) and characterized by the 
domains 

P(^;^)) = {/ + c « + e^^wi) I / G Co~(/), c G C} , 
where the functions wf{t) may be chosen as follows: 



wf{t)=wf\t) 



2.12) 



Remark. For a particular value of 7 the domain 'D{H'f^) is composed of functions, 
regular at t = 0. The corresponding SAE of Hi will be denoted by Hl'^^- 
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2.3 Resolvent 



The kernel Gk,i{t,t') of the resolvent of the radial Hamiltonian Hi satisfies the equation 

(^Mt) - k')G,,i{t, t') = ^^^^ - t') ■ (2.13) 
It basically means that if (^Hi — k^^ u = v for some u G V{Hi), then 

u{t) = I GkAt,t')v{t')dfit'. 



I 



In order to find the solution of the equation fl2.13p . consider the following ansatz: 
^ , J wl^^'^\t) wi^\t') for < t < f < 1, 



where the signs "+" and "— " should be chosen for / > and / < 0, correspondingly. It is 
clear that the function, defined by (12.141) . solves the equation fl2.13p for t t' and satisfies 
the boundary conditions of square integrability at the points t = and t = 1. (In the case 
|/| < 1 the requirement of square integrability at the boundary points is not sufficient to 
make the operator Hi self-adjoint; however, for such /, the ansatz (I2.14p also satisfies the 
regularity condition at t = and thus corresponds to the resolvent of the extension Hj'^^)- 
Taking into account the explicit form of the operator Hi, one may show that the 
required singular behaviour of the Green function at the point t = t' is guaranteed provided 
the condition 

dt'Gk,iit',t) 



t+o 1 

holds. Using ( 12.14^ . one can rewrite this condition as 



It follows from (I2.10p that the last relation is satisfied if we choose C^^ = C^J2. Substi- 
tuting this expression into (I2.14p . one finds a representation for the radial Green functions 
Gk,i(t, t')- 



3 Hamiltonian in the presence of a magnetic vortex 
3.1 Radial Hamiltonians 

Let us now add to the Hamiltonian the field of an Aharonov-Bohm magnetic fiux $ = 2711/, 
centered at 2; = 0: 



iv 


( dz 


dz 


2 


\ z 


z 



(3.1) 
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This choice of the flux position involves no loss of generality, since we have a well-known 
transitive SU{1, l)-action on D, which preserves the metric fl2.ip : 



Any gauge field configuration corresponding to a single vortex and a uniform magnetic 
field can be reduced to ^ = A^^^ + A^'"\ using the transformation (13. 2p combined with a 
gauge change. 

The Hamiltonian (12.21) in the presence of a vortex has thus the following form: 

r„ 1 „ 1 _ Mb ,^ . 46^ 



— w 



{drr + -dr + \id^ + ivf + {8^ + iv) - y^—^ ^4 • (3.3) 



This Hamiltonian still commutes with the angular momentum operator L. Radial Hamil- 
tonians H^^i are obtained by the restriction of Hy to the eigenspaces of L with fixed 
angular momenta / G Z. Namely, one obtains H^^i = Hi+v-, where the operators -ffagM are 
defined as in (12. 5p . Thus the only effect the AB vortex has on the formal Hamiltonians is 
the shift of the angular momentum variable by v. This observation allows to considerably 
simplify the derivation of many results, using the calculations from the previous section. 
Remark. As usual, for integer flux values some further simplifications occur. The Hamil- 
tonians iy(^) and are related by a gauge transformation 

= UH'-^^U^ U -.w^ e-^^^w;, (3.4) 

which is globally well-defined for G Z. The kernels of the resolvents of H^^^ and Hy in 
this case differ only by a factor of e'"^^'^~'^ \ and this change has no effect on the observable 
quantities. 



3.2 Self-adjointness 

From now on it will be assumed that — 1 < z/ < (it is clear from the above that this 
involves no loss of generality). Let us consider the full Hamiltonian Hy and restrict its 
domain to functions with compact support on the punctured disk: V^H^) = {D\{0}) . 
It was shown in the previous section that for |/| > 1 the operator Hi is essentially self- 
adjoint, and for |/| < 1 it has deficiency indices (1,1). One should then distinguish two 
cases: 

• u = . In this case has deficiency indices (1,1) and admits a one-parameter 
family of SAEs M^'^ with 7 G [0, 2n) and 

P(i^i^)) = {/ + c K + e^^^o ) I / G C^{D\{0}), c G C} . 

These Hamiltonians describe a purely contact (non-magnetic) interaction of a par- 
ticle with the AB solenoid. They have already been considered in [3] , so we will not 
pursue their study. 
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• — 1 < < . For such v the deficiency subspaces /C^ of the full Hamiltonian are 
generated by those of the operators and ii\j^u- Thus B.^ has deficiency indices 
(2, 2) and admits a four-parameter family of SAEs. Different extensions can be 
labeled by a unitary 2x2 matrix U and characterized by the domains 

I i=l,2 i=l,2 J 

where w^2 ortho normal elements of the bases of /C^, 

and II • II denotes the L^-norm on I with respect to the measure d^f 

Note that the diagonal matrix U describes magnetic point interactions acting separately 
in s-channel (/ = 0) and p-channel (/ = 1). Non-diagonal U introduces a coupling between 
the two modes so that the Hamiltonian no longer commutes with the angular momentum. 

Further analysis of spectral properties of H\l is a bit cumbersome in the general case 
(see, for example, the papers [21], [2], [TT], where such an analysis has been performed 
for the AB effect on the plane with and without magnetic field). We remark, however, 
that there exists a distinguished SAE of H^, whose domain consists of functions vanishing 
for t — > 0. This extension will be denoted by Hl^^. The next section is devoted to 

the calculation of its resolvent {^l^^ — k'^^ . The resolvent of any other SAE can be 
obtained from the latter using Krein's formula [1]. 

4 One-vortex resolvent 

4.1 Contour integral representations of the radial waves 

The main technical difficulty in the calculation of the resolvent kernel Gk{z,z') of the 
Hamiltonian Hl'^^ is the summation of radial contributions coming from different angular 
momenta: 

Gkiz, = E ^'^.'+'^(^' t'y'^''"''^- (4-1) 
Zir ^ — ' 

In order to address this problem, it is useful to introduce instead of the radial waves 
( I2.7p -( l2l9l) the functions depending on both t and ip: 

wl"(^) = (4.2) 

(I),,. _ r(x + b)T(x - b) ._a^+,i.„w 



W, [Z 



r(2x) 



<'{t), (4.3) 
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Combining these formulas with the relations (12.141) . (12. lip , one can rewrite the Green 
function (14. ip in the following way: 

—iu{ip—ip') , , 

GA^.z') = -^-^ (Ql:\z,z')+gt\z,z')) , (4.6) 

where the functions Q'j^\z, z') are given by 

glt\z,z') = Yl w«(.)w(^^'±)(/) for|^|>|^'|, (4.7) 

?(±)(2,/) = J2 ^l'''^\z)M'\z') foT\z\<\z'\. (4.8) 



The sums (14. 70 - (14. 8 p can be computed using a special set of solutions of stationary 
Schrodinger equation without AB flux, known as horocyclic waves |T2]. These solutions 
have the form 

1 - \z\ 



|2\X± 



^±{z,e) = — — _ (4.9) 

where 



1 / 1 

and 6 is an arbitrary complex parameter. Being considered as functions of 6, horo- 
cyclic waves \&±(2, 6') have an infinite number of branchpoints located at ^ = ±lnr + 
i{ip + IT + 27rZ). Let us introduce a system of branch cuts in the ^-plane as shown in 
the Fig. 1. The sheets of Riemann surfaces of the functions ^'^(z,^) are fixed by the 
requirement that the arguments of both 1 + z and 1 + ze^ are equal to zero on the 
line Im6' = ip. 

Recall that the Hamiltonians H^^^ and are related by the gauge transformation 
(13. 4p . Although this transformation is singular for non- integer values of the flux, one can 
still relate any solution of the equation {H^ — k^)w = to a solution of the same equation 
without AB field, (^(^) - k'^)^: = 0. However, since we have w = e ^'^'^ijj, the function ijj 
should be branched with the monodromy e'^'"^" at the point z = 0. Motivated by this well- 
known fact, we will try to represent radial wave functions (14. 2 1) - (14. 5 p as superpositions of 
elementary solutions (14. 9p . 

w(^) = / ^±{z,9)p{9)d9, 
Jc 

where C is an integration contour and p{6) is an appropriately chosen weight function. 
There will be three types of contours that will be important to us (see also Fig. 1): 
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Figure 1: Contours of integration in the 6'-plane 

• Contour C_|_(2;) starts at —oo + ia, surrounds the branch cut b+ = (— cxd + i(ip + 
tt), Inr + i{Lp + tt)] in a counter-clockwise manner, and goes to — oo + i{a + 2tt). 

• Contour C_(z) starts at cxd + i[a + 2tt), then goes counter-clockwise around the 
branch cut b_ = [— Inr + i{ip + vr), cxd + i(ip + vr)), and finally travels to oo + ia 
along the ray parallel to the real axis. 

• Contour Co (2;) joins two branchpoints: 61 = lnr + i{ip + TT) and 6*2 = — In r + i{(p + tt) . 

Real parameters a and 7 can be chosen arbitrarily; the only conditions they should satisfy 
are given by 

\ip — a\<7i, 0<7<— Inr. 

Assuming that Re/c^ < 0, one may now write a number of contour integral represen- 
tations for the radial waves fl4.2p - fl4.5l) : 



wj^) 






(4.10) 






JCo(z) 


(4.11) 


wj-^-) 




= ± [ ^+{z,9)e^U9, 

JC±{z) 


(4.12) 






= T [ '^+{z,9)e-^U9, 

Jc^{z) 


(4.13) 



where the functions '^±{z,9) are obtained from '^±{z,9) by replacing h — ^ —6. Although 
the validity of the representations ( I4.10p - (]4.13p can be checked directly, their general 
structure may also be a posteriori understood as follows. Consider, for instance, the 
functions W; (^) and w^'^'^'^'* (^) as defined by fl4.10p and f l4.12p . Continuation of these 
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functions along a counter-clockwise circuit enclosing the point z = amounts to simulta- 
neous shift of the branch cuts and integration contours upwards by 27r in the 6'-plane. This 
shift is in turn equivalent to simple multiplication of both functions by e^'^*'. Moreover, 
elementary solutions ( 14. 9p satisfy the relation 

L^±{z, 6) = {zd, - m)'^±{z, 6) = -de^±{z, 6), 

which means that RHSs of ( 14.10p and f l4.12p are common (multivalued) eigenfunctions of 
H^^^ and L, their angular momenta being equal to /. First function is regular for t — 1, 
since in this case the branch cuts pinch the imaginary axis. Similarly, the second function 
is regular for t — >■ 0. This implies (modulo constant factors that have to be found by a 
direct calculation) the relations (14.21) and (14. 4p . 



4.2 Summation 

Let us now turn to the calculation of the sums (14. 71) - (14. 8 1) . For simplicity the case \z\ > \z'\ 
is treated in detail and we only indicate the changes needed to handle another case. 
Substituting contour representations (14. 2p and ( 14. Sp into the relation ( 14. 7p . one obtains 

gl^\z,z') = T T I dej rf^2^-(^,^i)^+(^',^2)e'(^-^^). 

Since \z\ > \z'\, one may choose the contours C±{z') in such a way that 7^' > — Inr. 
Consequently, we have Re {9i - 62) < for all 9i E Co{z), 62 G C^{z') and Re (6^1 - ^2) > 
for all 6*1 G Co{z), O2 G Cj^{z'). Then it becomes possible to perform the summation 
inside the integrals and one finds 

gi-'\z,z') + gi-\z,z')= j de, j d92 vi/„(^,^^i)m/+(z',^^2) ^^^z^^^^ 

Co{z) C+{z')UC-{z') 

We would like to deform the contours C±{z') in the last integral over 62 so that their 
vertical parts compensate one another. Then C^{z') U C-{z') transforms into two hori- 
zontal lines, but one also earns a pole contribution coming from e^^ = e^^. Next, if we 
assume that ip — ip' ^ ivr, then the two lines can be deformed into Im6'2 = using 
quasiperiodicity in 62. Together with (14.61) . this leads to the following representation for 
the Green function: 

r ^-iu{v-^'+2n) q{0)^^^ ^ ^^(^^ if-ip'e (-27r, -tt), 

Gkiz,z') = I e-"'^^-^'^ Gf\z,z') + Ak{z,z') for ^ - ^' G {-n,n), (4.14) 
[ e-*^(^-^'-2-) Gf\z, z') + A,(z, z') for ^ - ^' e (vr, 2ti) , 

with 

Gf{z,z') = ^ I de ^^{z,e)i>^{z\e), (4.i5) 
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Co(2) Ime2=¥'' 

Similarly, assuming that \z\ < \z'\, one obtains an integral representation of the Green 
function which has exactly the same form as (I4.14p . except that the functions Gf\z,z') 
and Ak{z, z') are now given by 

Gf\z,z') = / de ^4z',e)^4z,e), (4.i7) 

Ak{z,z')= g^^^ y rf^i yrf02 ^-(^',^^i)*+(^,^^2)^^;z^73Y. (4.18) 

After some computations (technical details are outlined in the Appendix A) one may 
show that both representations coincide. Moreover, the integrals fl4.15p and fl4.17p can be 
carried out explicitly: 

Gf{z^z')={\:^}^\{u{z,z')), (4.19) 



where uiz, z') 



z' ~ z 



2 



has a simple relation with the geodesic distance between the 



l-zz' 

points z and z\ and the function C,{u) is given by 

c(.) = s^H^iiP^ (i - + "--^ - 1 - «) ■ (4-2°) 

Note that G^^\z,z') coincides with the well-known expression for the resolvent kernel 
of the Hamiltonian without AB field [9l [20]. This can also be seen directly from the 
representation fl4.14p . since Ak{z, z') in f l4.16p or (14. ISp obviously vanishes for u = 0. The 
function Ak{z, z') may also be written in a symmetric form: 

A,{z,z') = '-^ / de- J- ILLILI^] c(v{r,r',e)), (4.21) 

— oo 

with 

, ^, + r'^ + 2rr' cosh 6 , . . 

v{r, r', 6) = — . (4.22) 

1 + r^r'^ + 2rr' cosh 9 

In our opinion, the representation (14.140 and the formulas (I4.19p - (l4.22p constitute the 
most interesting results of the present paper. It is instructive to compare them with the 
known results in the flat space (cf. the relations (2.25)-(2.26) in [30] or the formula (5.10) 
from [27j). Notice that the 'free' part of the Green function is manifestly separated in 
(I4.14P from the vortex-dependent contribution Aj:{z,z'). 
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5 Spectrum and density of states 

The spectrum of the regular extension H^'^^ consists of three parts: 

• a continuous spectrum EE [(1 + '^h^)/ Fl? , oo); 

• a finite number of infinitely degenerate eigenvalues, which coincide with the usual 
Landau levels on the hyperbolic disk [9l |2l] in the absence of the AB field. These 
levels are explicitly given by 



i?2 



1 + 46^ -4 M6| - n- i 



(5.1) 



where n = 0, 1, . . . , rimax < \b\ — 1/2. Corresponding common eigenfunctions of the 
Hamiltonian H^^^ and the angular momentum operator L can be expressed in terms 
of Jacobi's polynomials (cf. the relation (13) in [M]): 

^Sl^.V') ~ tl'+H/2(i _t)|b|-"p^2|b|-2„-i,|i+H)(2t _ 1) e^'^. 
Here one should take / = 0, —1, —2, ... for 6 > and / = 1, 2, ... for 6 < 0. 

a finite number of bound states En'^ with finite degeneracy. The form of these 
eigenvalues depends on the sign of magnetic field. Namely, for 6 > one has 



i?2 



4b' -Alb -n 



(5.2) 



where n = 0, 1, 
be written as 



with n = 0, 1, 



71 

' max 



f>2 



< b — {u + 1) — 1/2. In the case 6 < 0, the eigenvalues may 

2" 



1 + 46^ -4 |6| -n + z/- - 



(5.3) 



n" 

' max 



< \b\ + u — 1/2. Common eigenstates of H^^^ and L are 



again given by Jacobi's polynomials: 



6 > : 
6 < : 



n,l 

{^,-) 
n,l 



(t, if) ~ t('+'^)/2(l _ i^b-n-{v+l) p{^b-2n~2{u+l)-l,l+v) - 1) e*'*^, 
(t,(^) ~ ^\l+u\/2^^_^^\b\-n+u p^2|feh2n+2z.-l,|Z+H)(2t_ 1) e^'^. 



For given radial quantum number n the allowed eigenvalues of the angular momen- 
tum are Z = 1, 2, . . . , n + 1 (for 6 > 0) and I = 0, —1, . . . , —n (for b < 0). 



Remark. The above expressions fl5.ip - fl5.3l) for the energy levels can also be extracted 
from the recent work [7]. It is worthwhile to emphasize that the discrete spectrum is 
absent for |6| < 1/2. 
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Let us now consider the density of states (DoS) on the hyperbohc disk. It can be 
obtained from the boundary values of the resolvent kernel on the real axis in the complex 
energy plane, using the following formula: 



p{E) = - ImTiGkiz, z' 

TT 



EG 



Both terms in the representation (14.141) of the Green function contribute to the DoS. 
The contribution of the free-resolvent kernel G''^ [z^ z') has been first calculated by Comtet 
[9j. His results (supplemented by an additional term [6j, coming from the discrete spec- 
trum) give the following expression for the DoS: 



P^°\E,z) = -lmGf\z,z' ^ z) 



k^=E+iO 



1 sinh27rA ^/^ 1 + 462" 



An cosh 2ttX + cos 2nb \ R? 



n=0 



Here 0fx) denotes Heaviside function and 



One can not expect that the DoS per unit area, induced by the AB field, will also be 
constant on D. However, it should depend only on the geodesic distance between a given 
point on the disk and the flux position. Indeed, since the function A(z, 2;') is non-singular 
for z — > z', the vortex-dependent part of the DoS is given by 

p('^)(E,^) = -ImAfc(t) , (5.5) 

where the function Afc(t) is obtained from A(z, 2;') by setting = = r'^ = t: 

, sinvrz/ f ,^ e(i+^)'' fl + te-^\\f 2t(l + cosh^) \ 

—00 

As it stands, the representation (15. 6p is valid in the left half-plane Re k"^ < 0, where the 
function Ak{t) is analytic. However, the DoS is determined by the singularities of Ak{t) 
that occur on the positive part of the real axis (we may expect there a finite number of 
poles and the branch cut [(1 + 4iP')/ E? , 00), corresponding to the continuous part of the 

spectrum of Hl^^). One could try to construct the appropriate analytic continuation of 
Afc(t), considering (15. 6p as a contour integral and then suitably deforming the contour. It 
seems, however, that this approach does not lead to any satisfactory result because of the 
complicated singularity structure of the function under the integral sign in the 6'-plane. 
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An alternative method consists in the following. Remark that the vortex-dependent 
contribution to the DoS in the whole hyperbolic space 



P 



{E)= [ diip^^\E,z) 

J D 



(5.7) 



has a finite value, since 



1 +te 
1 +te 
1 

An 



C 



2t(l + cosh I 



1 +t2 + 2tcoshe^ 

ln2t + ln(l + cosh6') + 27i5 + ?/'(x + &)+^(x-&) + o{l) for t ^ 0, 



1 T{x + h)V{x-h) 



, 47r 



r(2x) 



'l + e^)x+b(l + e-e)x 



— + o ((1 - t)'^) for t ^ 1. 



If one now integrates Afc(t) over spatial coordinates (see Appendix B) and then consid- 
ers the analytic continuation of the result to the complex energy plane, the following 
expression for p^'^\E) can be obtained: 



p(")(E) = -— Im 



(X - 6 + z/) i){x -b) - ipix -b + u + 1) 



{5.t 



4vr 2x - 1 
+ {x + b-i^~l)\Hx + b) ~ij{x + b-i^-l)\ , , 

L J J k^=E+iO 

where the contributions of the discrete and continuous part of the spectrum are given by 



£(n + l)5(E-E^+))- £(n-z.)5(E-E(°)) 



for 6 > 0, 



n=0 
n" 

max 



n=0 



J] (n + 1) 5 - El^"'-^) - ^ (n + + 1) 5 - E^^^) for b < 0, 

(5.9) 



^ n=0 



n=0 



pM{E) = -^i-Q{E~ 



1 + 452 \ ( X sinh 27rA + {\-b + v) sin 27r(6 



cosh 2tt\ + cos 27r(6 — v) 
X sinh 2-nX + {\ — b + u) sin 27r6 



(5.10) 



cosh 27rA -|- cos 27r6 

and the parameter A is defined as in (15. 4p . 

At last we add a comment concerning the fiat space limit {R — >■ 00) at zero magnetic 
field (6 = 0). In this case the representation (15.81) for the vortex-dependent DoS transforms 
into 



(5.11) 
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This result has been first obtained in [10], and it has important consequences in the theory 
of disordered magnetic systems (see, for example, [T3|, [HI [I5]). Obtaining the relation 
(15. lip directly from (15.101) is more subtle; one should consider p^\E) as a distribution 
and supply it with a proper regularization at the edge of the spectrum, i. e. as A — > +0. 

6 Discussion 

We have studied the Hamiltonian of a particle moving on the hyperbolic disk in the 
background of a uniform magnetic field and the AB gauge potential. The density of 
states and the resolvent of this operator have been calculated in a closed form, using 
Sommerfeld-type integral representations for the radial waves. 

The above discusssion does not exhaust all problems related to the system under 
consideration. First of all, there remains a technical question of analytic continuation 
of the representation (14.211) to energy values with Re fc^ > 0. Such continuation would 
allow to investigate the curvature dependence of various vacuum quantum numbers, e. g. 
fractional charge, magnetic flux and angular momentum (see, for instance, [33]), induced 
by the AB vortex. One may also try to address the latter problem, applying the technique 
we have used in the calculation of the density of states. 

Another interesting question concerns the generalization of our results to the case of 
the Dirac Hamiltonian. This problem, in its turn, appears to be non-trivially related 
to the theory of isomonodromic deformations on the hyperbolic disk [281 [311 [32] . More 
precisely, it was shown in [31] that the two-point isomonodromic tau function of the Dirac 
operator on the Poincare disk provides a solution to Painleve VI equation. The explicit 
form of this Painleve transcendent in a particular case of zero magnetic field has been 
conjectured by Doyon [17j, whose idea was to replace the tau function by its physical 
analog — a correlator of monodromy fields, and to sum up the corresponding form factor 
expansion. The missing ingredient in the rigorous proof of this result and its generalization 
to the case of non-zero field is the Green function of the Dirac operator with one branch 
point, precisely analogous to the resolvent found in the present paper. Moreover, the 
knowledge of this Green function enables one to compute all form factors of monodromy 
fields on the hyperbolic disk [T7], in particular, those of Ising spin and disorder fields 
[T5[ We leave these problems to a future publication. 

Acknowledgements. I warmly thank Dublin Institute for Advanced Studies, where 
a large part of this work has been done, for hospitality. I am very grateful to Sasha 
Povolotsky for useful conversations and help with the literature. 

Appendix A 

Here we describe a method of obtaining the expressions (I4.19p - (l4.22p for the Green func- 
tion from the integral representations (14. 150 - 04.180 . First, note that horocyclic waves 
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fl4.9p satisfy the following identity: 



e) = {a + (3e'') [a + (3e') "^^"^ ' ^^{s, 0, (A.l) 

where we have introduced the notation 

and assumed that |Q;|^ + |/3p=l. Let us now make the change of variables flA.2p in (14. 151) 

az' - (3 

(introducing similarly the variable s' = — = instead of z'). It is then straightforward 

—f3z' + a 

to check that fl4.15p reduces to 

Gr(-') = ^(|ii^)7|;i^)7 d(^-(M)'KW,0 (A.3) 

47r\/3s + aJ \(3s' + aJ Jc^is) 

We may now adjust the parameters a and /3 in such a way that s' = 0. For example, one 
can take 

a = — -^=^= , p 



\zV ' ' V^- 



z — z' 



Since in this case s = —, the relation (1A.3P transforms into 

1 — zz' 



Co(s) 

I f TK- (A.4) 



1 fl-Zz'\ , 2\X- /" 



47r VI - zz' ) ' ' ' ' ^ (1 - |s|e-€) 

ln|s| 



After the change of variable = ^ in the last integral, we obtain standard 

\s\ 

representation of Gauss hypergeometric function. 



(lei) (-NT/ 



'''' ^Tl\\--ZZ')^ ' ' ^ 7 (l_(l_|s|2)i)X+;' ' 



which immediately gives fl4.19p - fl4.20p . Performing analogous manipulations with the 
representation f ]4.17p . one finds the same answer. 

Let us now consider the representation (14. 160 for the function Afc(z, z'\ Interchanging 
the order of integration and then introducing instead of a new variable 62 = O2 ~ Oi, 
one can check that this relation transforms into 

A,{z, z') = -^r— e-^-^^-'^ ) / dh F{z, z' , h) , (A.5) 

z-Ki J - 1 

Im 62=^' —i^—T^ 
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where the function F{z, z', 62) after additional change of integration variable 61 ^ 61 = 
61 + 62 can be written as 

F{z, z', 62) = ^ / d9i ^^{z, h - ~e2)^+{z, h) . (A.6) 

In the last expression, Cq^{z) represents the integration contour obtained from Cq{z) by 
shifting it by 6*2. 

One may now use the trick described above to eliminate the function from 
flA.6p . Namely, consider further change of variables: 



z'e^i + 1 



Somewhat tedious but fairly routine calculation shows that the integral ( ]A.6p . being 
rewritten in terms of ^, reduces to 

F{z. z\ O2) = j- (\±4^) ' (1 - ssr [ , _1 , , (A.7) 



where we have introduced the notation 



r e 



Re 62 _j_ ^ ^g— Re92 _j_ ^' 



e*^, 5 = ^e"*'^. (A.8) 

][ _|_ g Re 62 ]^ _j_ g— Re 62 

Contour denotes the horizontal line segment in the ^-plane, joining the points In |s| +i<y9' 
and — In |s| +iip'. Rewriting flA.7p in terms of ordinary integrals and using flA.Sp . one gets 
a relation analogous to (]A.4p : 

1 / 1 _|_ rr' p Re 62 \ 

F{z,z',e2) = \ s (i-sSr-x 

47r \ 1 + rr'e-^'^^ 



X 



-iln|ss| 



]^ _j_ j^^' g Re 02 
1 _|_ ff' g~I^e §2 



C{\s~s\). (A.9) 



Finally, combining the last formula with ( lA.SI) and ( lA.SI) . we find the relations ( ]4.2ip - 
f l4.22p . The same result can be obtained from the representation (14.181) in a completely 
analogous manner. 
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Appendix B 

In order to integrate the function Afc(t) over the hyperbohc disk, consider the following 
identity: 

b 



Vl + te^ J ^ Vl + t2 + 2tcosh^y 

oo 

(1 - tf^ f de' 



47r I [(1 + e-e)(l + e-^') - (1 -t)e-e]^"^ 
1 



X 



[(1 + e^)(l + e^') - (1 - t)(e^ + e^' + e^+^')] 



Substituting this relation into the integral / d/i Afc(t), integrating over the angle 9? and 

J D 

introducing instead of t a new variable s = 1 — t, one obtains 



\ rf/i Afe(t) = ^ / de / dff ^—^ Tie, (B.i) 



[(1 + e-^)(l + e-^') - se-^]^-' [(1 + e^)(l + e^') - s(e^ + e^' + e^+^')]''^' 



5 

After further change of variable s ^ u= — 57— the function Tie, e') 

s + (1 + e'^)(l + e*^ )(1 - s) 

can be rewritten as 



T{e,e') 



u 



^^-^du 



(1 + e<^){l + e^') 7o [e-^-^' + (1 + e-'')uf- 



Integrating once by parts, one finds 



{2x - 1)(1 + e^) Jo [e-'-'' + (1 + e-^>] 



x-6+1 



(B.3) 



Let us consider the contribution of the term flB.2p to the integral fIB.ll) . Introducing 
instead of e' a new variable f = [1 + + e~^~^ , we may write 

poo poo r,{i+y)0 1 /•! poo p{i+y)e ?>x-^'-i 

/ de de' -TAe,e') = / dv de — -. (b.4) 

In the contribution of the term ( IB. 31) . we first perform the integration over then replace 
e — >■ —e and finally exchange the order of integration over e and u. The result looks quite 
similar to (]B.4p : 

poo p(l+u)e 1 pi poo -ue yX+b-l 

de de' - — -^2(^,^0 = 7^ — T du de- — — -. (b.s) 



00 



l + e^ ' ' 2x-Uo J-00 (l + e^yi + ue- 
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The integrals over 9 in (]B.4p - (lB.5p can be easily calculated by residues. For example, one 
obtains 



Sm TTZ/ 



de 



v^^" - I + {I + u){l - v) 



Subsequent integration over v and u leads to the final result: 



D 



dfiAkit) = --^——^^{x-b + iy)^^{x-b)-^{x-b + u + l) 
+ {x + b-i^-l) k(x + b)-ij{x + b-l-u) 



+ 



where iJj{x) denotes the digamma function. 
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